Jamming and asymptotic behavior in competitive random parking of bidisperse cars 
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We propose a generalized car parking problem where either a car of size a or of size ma (m > 1) 
is sequentially parked on a line with probability q and (1 — q), respectively. The free parameter 
q interpolates between the classical car parking problem at either extreme (q — and q = 1) and 
the competitive random sequential adsorption of a binary mixture in between. We find that the 
coverage in the jamming limit for a mixture always exceeds the value obtained for the uni-sizcd 
case. The introduction of a bidisperse mixture results in the slow approach ( ~ t -1 ) to the jamming 
limit by the smaller species while the larger species reach their asymptotic values exponentially fast 
~ t~ 1 e~ i - m ~ 1)qt . 

PACS numbers: PACS number(s): 05.20.Dd,02.50.-r,05.40-y 



The adsorption of particles on a flat substrate is a com- 
mon phenomenon and it is of wide interest in physics, 
chemistry, biology and in many other branches of science 
and technology. Examples are the adsorption of macro- 
molecules .colloidal particles, bacteria, proteins and latex 
particles |l|, ||, Hi- Due to its wide range of applications, 
there has been continuing effort to develop and enrich 
our understanding of the topic through all the avenues 
of research comprising experimental, numerical and ana- 
lytical means GJ. 

From the theoretical viewpoint, the simplest model 
that can still capture the generic features of the adsorp- 
tion phenomenon is the kinetics of random sequential 
adsorption (RSA) of monodisperse particles on a clean 
d-dimensional substrate In this process, particles 
are deposited randomly, one at each time step, with the 
restriction that overlap is forbidden. The one dimen- 
sional continuum version of the RSA process is popu- 
larly known as the random car parking (RCP) problem. 
This simple RSA model ignores diffusional relaxation, 
precise particle-particle and particle-substrate interac- 
tion (though the hard-core interaction is incorporated) 
and desorption or evaporation of the adsorbed particles. 
Nevertheless, there has been a continued effort to bridge 
the gap between the experimental conditions and the sim- 
plest form of RSA. As a result, there are many interest- 
ing variants of RSA. Examples include the cooperative 
sequential adsorption ||, the accelerated random se- 
quential adsorption 0, the ballistic deposition ||, the 
RSA on disordered substrates the RSA of growing 
objects (T(|. Several attempts have been made to in- 
clude the transport of deposited particles by diffusion 
(see and references therein). One issue inherent to 
the experimental realization is that the adsorbing parti- 
cles are not always perfectly monodisperse. Polydisper- 
sity is almost an inevitable property in many real life 
situations. The binary mixture is the simplest and the 
first step toward the understanding of polydisperse sys- 
tem. There has been an increasing interest in the study 
of the RSA of mixtures with a strict lower and upper size 
cut-off jl2|, [H| [li], , and continuous mixtures obeying 



a power-law size distribution fllq | . It is important to note 
that in the former case the structure of the monolayer in 
the long time limit is described by the jamming cover- 
age, whereas in the latter case the resulting monolayer 
is uniquely characterized by the fractal dimension of the 
arising pattern jl6|. Much of the theoretical effort was 
restricted only to a binary mixture of particles with very 
large size differences However, mixtures with a 

narrow distribution of sizes have more practical relevance 
than mixtures with large size differences. 

In this article we consider the RSA of two species and 
incorporate a parameter that controls the competition 
between the adsorbing species. Hence we shall call it 
competitive random sequential adsorption (CRSA) [ fl7j . 
The motivation of setting the present problem came from 
its potential applications in polymer, aerosol and col- 
loidal science. Consider the case of the adsorption of 
polymer chains onto a surface or the reactions along poly- 
mer chains, which was in fact the motivation of the orig- 
inal work by Flory Q . The polymer chains are naturally 
not all of the same length, instead they are a mixture 
of chains of lengths differing by an integer multiple of 
a certain size unit. Of course, the ideal case to consider 
would be the study of a polydisperse mixture with strictly 
lower and upper size cut-off which would warrant greater 
complexity. Nevertheless, the present work still captures 
all the essential and generic aspects of the adsorption of 
mixtures in the simplest possible way. Thus, our model 
provides useful insight for the analysis of more complex 
problems. We find that the mixtures exhibit a behavior 
which is qualitatively and quantitatively different from 
that of the monodisperse adsorption. 

It is worth to mention that although the RSA prob- 
lem is based on a simplified picture, yet it is analytically 
tractable only in one dimension. This is typical for many 
problems of statistical physics. The underlying reason is 
that the shielding property can only be found in Id. Here 
shielding means that in Id the adsorption of a particle 
on an empty interval separates it into two disconnected 
smaller intervals, both having the same geometry as their 
parent interval and hence they can be treated as an in- 
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dependent interval. Therefore, much of the information 
in higher dimensions is provided by numerical simula- 
tions or by experiment. We can identify two separate 
characteristic time scales in adsorption phenomena: the 
time between depositions, td, and the time required for 
the deposited particles to reorganize themselves on the 
substrate, t r . In the context of the present work, we 
restrict ourselves to the case where the relaxation time 
is very large, td « t r , so that the deposited particles 
on the substrate do not have enough time to reorganize 
to find the minimum energy configuration and hence the 
system is driven out of equilibrium. Such a situation is 
expected when the attractive contact energy between the 
adsorbed particles and the substrate is sufficiently higher 
than the thermal agitation energy ~ fcsT. Restriction 
to the monolayer formation and the irreversibility of the 
process make it one of the simplest yet interesting and 
challenging problems in statistical physics. Despite the 
inherent simplicity of the RSA model, it has in fact re- 
produced many experimental results. Namely, the RSA 
of monodisperse hard particles was first shown to de- 
scribe the experimental results of the adsorption of pro- 
teins [|l9|| , and later that of latex and submicron colloid 
particles [^0| on flat surfaces, for example. In these pro- 
cesses, it is known that the relaxation is typically very 
slow and hence the simple RSA algorithm works well. 
The known key features of the RSA processes are: 

a) The process results in a unique jamming coverage 
in the thermodynamic limit. The final monolayer 
is uniquely characterized by the jamming coverage 
6d(oo) in d dimensions, defined as the fraction of 
the substrate covered by the depositing particles, 
which is always less than unity due to its irre- 
versible character, and the restriction to the mono- 
layer growth. 

b) Near the jamming limit, the approach of the cover- 
age 9d(t) to its final value 9d(oo) obeys the relation 



Uoo)-e d (t) -i- 1 '*. 



(i) 



This is widely known as Feder's law ||l9[ and the 
theoretical arguments supporting it have been pre- 
sented independently by Swendsen and Pomeau 
@. 

c) Another surprising result in RSA of d-dimcnsional 
hyperspheres is that an estimate of the final jam- 
ming coverage 9d(oo) can be obtained with a rea- 
sonably good degree of accuracy from the result of 
the one dimensional model (RCP) through the re- 
lation 



9 d (oo) « (#iM) d . 



(2) 



Although the strict equality as conjectured by 
Palasti pH is not valid, the relation often provides 
a remarkably accurate estimate. 



The problem of the competitive random sequential ad- 
sorption can be formulated as follows. Assume that we 
have a particle reservoir which contains a binary mix- 
ture of chemically identical particles of two different sizes. 
Further, assume that a fraction q of the particles in the 
reservoir is of size a and the rest, i.e. a fraction (1 — q) 
of size mer, m > 1. The algorithm of the process is then 
as follows 

i) A particle is picked randomly for deposition on the 
substrate: it can be of size a with probability q or 
of size ma with probability p = 1 — q. 

ii) Randomly a position is chosen on the substrate. 

Hi) The particle is deposited on the substrate if the 
chosen position falls on an empty interval and it is 
large enough to accommodate the particle. Other- 
wise the particle is rejected. 

iv) In either case the time is increased by one unit. 

v) The process is repeated until no more gaps with 
size larger than a are left on the substrate. 

The task now rests on translating the problem in ques- 
tion into a set of equations for the gap size distribution 
function G(x, i) to describe how it evolves in time. Before 
doing so, we first note that the two events (adsorption of 
a and ma) are mutually exclusive, i.e. at each time step 
only one particle is adsorbed. This is ensured by the 
fact that either the particle of size a is adsorbed with 
probability q or a particle of size ma is adsorbed with 
probability (1 — q). Then, we can immediately generalize 
the classical car parking problem and write the following 
set of equations for G(x, t) for x > ma 



dG(x,t) 

m 



)G(x,t)+2q G(y,t)dy 

J x-\-cr 

pOO 

2p / G(y,t)dy, (3) 

J x-\-ni(T 



and for x < ma 

dG{x,t) 
dt 



poo poo 

2q G(y,t)dy + 2p G(y,t)dy 

J x-\-cr J x-\-7na 



H(x — a)q(x — ct)G(x : t) : 



(4) 



where a av = {m + q(l — m)}a is the average particle size 
and H(z) is the Heaviside step function. Equations (|^) 
and (||) satisfy the conservation law 



(x + {m + q(l - m)}a^jG(x,t)dx = 1. (5) 



To solve the kinetic equation, we seek in Eq. ^ a trial 
solution of the form 



G(x,t) = i4(t)e-<*- a "> B W 



(6) 
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where A(t) and B{t) are yet undetermined functions. 
Substituting the trial solution into Eq. (||), we obtain 
the following two differential equations for A(t) and B(t) 



din Ait) e" 



2p- 



dB(t) 



= 1. 



(7) 



dt t t dt 

Solving these equations subject to the initial conditions 



lim 

L >OG 



\ G(x, 0)dx = 0, lim / xG{x,t)dx = l, (8) 
Jo 1 >°Jo 

and satisfying the conservation law Eq. (Q), gives 



A(t) = t 2 F(at) and B(t) = t, 



where the dimensionless quantity F(at) is 



F(at) 



= p -2./ CTt {-3(l-e — )+p(l-e- 



1 )}du/u 



(9) 



(10) 



Note that if we set a = in Eq. (Q) , we can still recover 
the well known solution of random scission model of a 
binary fragmentation process G(x,t) = i 2 exp[— xt] [p3| . 
The complete solution to Eq. (||) is therefore 



G(x,t) = t 2 F(at)e- (x - a ^ )t for x > ma. 



(11) 



We can solve Eq. ([|) directly for a narrow size difference 
(1 < to < 2), since in this case only the known solution 
( |TT| ) contributes to the integrals in (^). To this end we 
seek a trial solution of the form 



G(x,t) = C{x,t)e- q( - x -^ 1 . 



(121 



In this domain, a is the minimum gap size, and hence the 
minimum value of the lower limit of the first integral of 
Eq. (Q) is 2a. Therefore, if we restrict to 1 < to < 2 we 
can substitute the solution of Eq. (|j) into both integrals 
and obtain the solution for a < x < ma 



G{x,t) 



F(as)s(2qe 



r(mp—l)s 



+ 2pe~ 



x e -{x{qt+ps)-qat)^ s ^ 



(13) 



Defining G(x, t) = G(x — a av , t) for x > ma we find that 
the solution of Eq. (^) satisfies the identity 



G({3 



,}\,t/\) = \- 2 G(x-a av ,t). (14) 



This relation is the hallmark for the existence of scale in- 
variance and it is equivalent to the data-collapse formal- 
ism p4| . It means that if the deposition rate (x — a av ) is 
increased by a factor A, and the observation time is de- 
creased by the same factor, then the resulting structure 
would look the same except for the numerical pre-factor. 

From the gap size distribution function one can im- 
mediately find the total coverage #t(£) by the adsorbed 
particles through the definition 



it) 



i 



x G(x 7 t)dx. 



(15) 



For < q < 1, both species contribute to the total cover- 
age 8r{t) = 6s (t) + In order to obtain expressions 
for the coverages by the small and the large particles, 
6s(t) and #i(i), we substitute the solutions of Eq. (||) 
and (Q) in 0T{t)- After some lengthy yet simple manipu- 
lations, we obtain the contribution of the small particles 
0s(t) 

9 s (t) = 2 f dsF(s)se kps (qe- s +pe- {k+1)s ) 
Jo 



-k(ps+qt) 



1 



-ks 



1 



(ps + qt) 



+q / F(s)(l + ks)e~ kqs ds, 



(16) 

where k — (to — 1), and the contribution of the large 
particles #z,(t) 

9 L (t)=pm[ F{s)e- kqs ds. (17) 
We calculate the coverage in the jamming limit (t — ► 




FIG. 1: Coverage in the jamming limit vs q for the respec- 
tive contributions from small and large particles for the case 
m = 1.5. Symbols repre sent simulational results, lines the 
analytical solution (P_6Ul7j) . Simulations were carried out for 
a substrate length of 10000 small cars, each symbol is an av- 
erage over 20 independent realizations. 

oo) of our model by direct numerical integration of the 
above two equations. The two terms in 9s are the con- 
tributions made by the adsorption of small particles in 
the domains a < x < ma and x > ma, respectively. 
The individual concentration of cars of size a and ma on 
the monolayer depends on their respective probabilities 
q and p. The final concentration of small particles will 
always exceed q, since the smaller particles have a higher 
success rate in the competition with the larger species, 
in particular when the substrate is already densely pop- 
ulated. For to = 1, the first term in 9s{t) is zero and the 
remaining terms of 9t (t) yield 



9 T (t) 



F{s)ds = R (t), 



(18) 



where #R,(i) is the coverage for RCP model with 6*r(oo) 
0.748..., the famous Renyi- limit. 
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To test our analytic results, we simulate the CRSA 
problem on a computer. Naturally, the simulations re- 
strict on finite substrate lengths L. However, for suffi- 
ciently large lengths L > a the effects of finite L are 
small. The algorithm we use follows literally the steps 
i) - v) described in the introduction. The time scale in 
the simulations is set by the number of particles brought 
to the substrate (if adsorbed or not). The discrete 'loop 
index' i of the simulation relates to the time scale of the 
rate equations t as 

t = i^. (19) 




q 



FIG. 2: The total coverage vs q for various values of m. Di- 
amond symbols refer to simulations with a substrate length 
of 10000 small cars, and for the case m = 1.95 results from a 
simulation with L — 5000 are shown for comparison as square 
symbols.. 

Using this time scale and expressing all lengths in the 
simulation in terms of a, the simulational results can be 
directly compared to the solution to the rate equations. 
In Fig. [l| the results for the jamming coverage are shown 
for the case m = 1.5, where the parameter q is varied. In 
general, we find an excellent agreement between simula- 
tion and our theory. The total coverage is found to be 
always larger for binary mixtures (q 0, q 1) than for 
the uni-sized case (q = or q = 1). For a given m the 
maximum coverage #™ ax (oo) occurs when the concentra- 
tion of the smaller species in the reservoir is small, i.e. at 
q — ► + . In this limit, the larger species are adsorbed 
first until their jamming limit is reached. Then the re- 
maining gaps are filled by the smaller species. In other 
words, the two processes (adsorption of small and large 
particles) are decoupled. As a result, the chronology of 
the adsorption of the small and large species strongly 
influences the dynamics of the system vis-a-vis the fi- 
nal coverage as well as the final number density of each 
species on the substrate. Note that the particles are 
drawn from an infinite reservoir, i.e. no matter how small 
the concentration of small particles in the reservoir will 
be, they still may be picked up at the expense of many 
failures. The maximum coverage for a given m occurs at 
q — ► + and it decreases monotonically as m decreases, 
eventually coinciding with the Renyi- limit at m = 1. 
The highest coverage is obtained for a large size ratio 



(m — > oo) of the binary mixture (see Fig. g). In the 
limit m — ► oo and a > the maximum total coverage 
is 

6^ ax (oc) = R (oo) + (1 - r (oo))0r(°°) = 0.937. . . . 

(20) 

The simulation with m = 20 shown in Fig. |2| is already 
very close to this limit. This implies that a maximum 
of 93.7 percent of the total space can be covered by a 
randomly adsorbed binary mixture. Note that the cov- 
erage for the CRSA of a binary mixture has no smooth 
transition to that of a mixture of points and finite sized 
particles in the limit m — ► oo; the small particles, no 
matter how small their finite length is, will always con- 
tribute to 9s, and thus, the total coverage 9t is always 
larger than the Renyi-limit (9r(oo). Only if their length 
is indeed zero, so that 6*5(00) = 0, the only role the 
point particles can play is to prevent the large particles 
from reaching the Renyi-limit. Also this problem may be 
solved exactly J25| . The result for the total coverage 

9 T {t)=p f F{s)e~ qs ds. (21) 
Jo 

In this case the total coverage stems from the large par- 
ticles alone and Feder's law does no longer holds. 




0.001 0.010 0.100 1.000 10.000 100.000 1000.000 

TIME 

FIG. 3: Time development of the coverage for q = 0.5, 
m = 1.5. Lines (solid, dashed, dotted) are theoretical results, 
symbols from simulations. 

From the gap size distribution function one can also 
obtain the number density of the individual species, using 
the definition 

poo 

n s + n L = / dyG(y,t). (22) 
Jo 

Using similar manipulations as we did for the jamming 
coverage we find 

ns(t) = 9 s (t) and n L {t) = 9 L {t)/m. (23) 

Fig. H shows the time dependence of the coverage by 
small, large, and by both species. 0i,{t) approaches its 
asymptotic value (the jamming limit) much more rapidly 
than 9s{t). In order to quantify the approach to the 
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respective asymptotic values we plot ln[0,s(oo) — #s(t)] 
vs hi[t] (Fig. |j). The small particles reach their jam- 
ming limit algebraically, obeying Feeler's law (Eq. (El)). 
Asymptotically the coverage for the large particles, Eq. 
( |l7| ) behaves as 

0£(t)~0 L (oo)-e-< m - 1 >«7t. (24) 

Thus the large particles reach their asymptotic value ex- 
ponentially with an algebraic prefactor hence not 
obeying Feder's law. This implies that the asymptotic 
behavior of the total coverage is dominated by the slow 
dynamics of the smaller specie, approaching the final val- 
ues algebraically as t~ x (see Fig. [|). Thus, we verified the 
validity of Feder's law for the present system. Note that 
according to Eq. ( p3|) the number densities will follow 
the asymptotic behavior of the respective coverages. 

The present work provides an exact theoretical basis 
for the RSA of mixtures. Our results are in good qualita- 
tive agreement with the numerical simulations presented 
by Meakin and Jullien |l2[ for the two-dimensional CRSA 
of disks. Especially, they find a more efficient coverage 
for the mixture as well. This is in contrast to the study 
of CRSA by Bonnier Il3] in Id, who claims that this 
behavior was typical for lattice models, while his con- 
tinuum model predicts a less efficient coverage by the 
mixture than by the RSA of uni-sized objects. The con- 
tinuum model we have studied in this work is consistent 
with its lattice counterpart. Our results are supported 
by the good agreement with the direct simulation. Fur- 
thermore, the dependence of our results on q and to is in 
excellent qualitative agreement with the results reported 
in yl|. Further, Meakin and Jullien too confirm Feder's 
law for the small particle species, and find an exponen- 
tial approach of the coverage by the large ones. In this 
context it is interesting to verify Palasti's conjecture (Eq. 
(0)), taking the values from Fig. 3 of the paper by Meakin 
and Jullien for the jamming coverage by a binary mix- 
ture of disks. We find a deviation of at most 5% from 
the estimates provided by our one-dimensional data and 
Eq. (||). This confirms that Id models, despite their sim- 
plicity, indeed may prove useful in gaining information 
about higher dimensions, in contrast to many problems 
in equilibrium statistical physics. There, in general, Id 
systems often have very little relevance to 2d or to higher 
dimension, e.g. no thermodynamic phase transitions are 
possible in Id. In RSA there is no phase transition in 
the adsorbed monolayer regardless of its dimensionality. 
The Id model already has all the ingredients such as to- 
tal irreversibility with the associated memory effect or 



exclusion or blocking phenomena that make their behav- 
ior nontrivial and qualitatively similar to that observed 
in higher dimensions, except a trivial quantitative dif- 
ference. Of course, the RSA with hard objects covers 
less and less fraction of the substrate as the dimension 
increases. 

In conclusion, we have shown that an analytical model 
for the CRSA of a binary mixture in one dimension pre- 
dicts a jamming limit that is always larger than that for 
the uni-sized case. In the limits q = 0, q = 1, as well 
as rn = 1 we recover the classical result 0r = 0.748 for 
the equal sized RSA process. In all other cases, i.e. for 
true binary mixtures, we find a larger jamming limit in 
theory and simulation. We have explicitly shown that 
not only the size ratio between the two species to which 
is an important parameter, but also the concentration 
of the individual species in the reservoir which depends 
on q is an equally important parameter. It appears that 
the smaller the concentration of the small particles in the 
reservoir, the higher the jamming coverage. In general, 
our analytical results are in excellent agreement with di- 
rect numerical simulations. In addition, they show an 
excellent qualitative agreement with numerical simula- 
tion results in higher dimension which emphasize the im- 
portance of Id model. The time-asymptotic approach 
of the jamming limit is dominated by the contribution 
of the small particles, and we confirm the 1/t behavior 
predicted by Feder's law in one dimension. The large 
particles reach their contribution to the jamming limit 
exponentially with an algebraic pre-factor. 




TIME 

FIG. 4: The approach of 0s(t), and thus Orit), to its asymp- 
totic value obeys Feder's law (equation ([j])). 
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